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Atomic and molecular adsorption on transition-metal carbide (111) surfaces from 
density-functional theory: A trend study of surface electronic factors 
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Chalmers University of Technology, SE-^12 96 Goteborg, SwedeT\j 
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Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark 

This study explores atomic and molecular adsorption on a number of early transition-metal car- 
bides (TMC's) by means of density-functional theory calculations. The investigated substrates are 
the TM-terminated TMC(lll) surfaces, of interest because of the presence of different types of sur- 
face resonances (SR's) on them and because of their technological importance in growth processes. 
Also, TM compounds have shown potential in catalysis applications. Trend studies are conducted 
with respect to both period and group in the periodic table, choosing the substrates ScC, TiC, 
VC, ZrC, NbC, <5-MoC, TaC, and WC (in NaCl structure) and the adsorbates H, B, C, N, O, F, 
NH, NH2, and NH3. Trends in adsorption strength are explained in terms of surface electronic 
factors, by correlating the calculated adsorption energy values with the calculated surface electronic 
structures. The results are rationalized with use of a concerted-coupling model (CCM), which has 
previously been applied succesfully to the description of adsorption on TiC(lll) and TiN(lll) sur- 
faces [Solid State Commun. 141, 48 (2007)]. First, the clean TMC(lll) surfaces are characterized 
by calculating surface energies, surface relaxations, Bader charges, and surface-localized densities 
of states (DOS's). Detailed comparisons between surface and bulk DOS's reveal the existence of 
transition-metal localized SR's (TMSR's) in the pseudogap and of several C-localized SR's (CSR's) 
in the upper valence band on all considered TMC(lll) surfaces. The spatial extent and the dangling 
bond nature of these SR's are supported by real-space analyses of the calculated Kohn-Sham wave 
functions. Then, atomic and molecular adsorption energies, geometries, and charge transfers are 
presented. An analysis of the adsorbate- induced changes in surface DOS's reveals a presence of both 
adsorbate-TMSR and adsorbate-CSR's interactions, of varying strengths depending on the surface 
and the adsorbate. These variations are correlated to the variations in adsorption energies. The 
results are used to generalize the content and applications of the previously proposed CCM to this 
larger class of substrates and adsorbates. Implications for other classes of materials, for catalysis, 
and for other surface processes are discussed. 

PACS numbers: 68.43.Bc, 73.20.At, 73.20.-r 
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I. INTRODUCTION 

Many studies of TMX's (TM = transition metal and X 
= C or N) are motivated by a curiosity on the properties 
of TMX's, such as their mixture of covalency, ionicity, 
and metallicity, and by a suggested importance for het- 
erogeneous catalysis. The attention to TMX's as poten- 
tial catalysts started with an observation by Levy and 
Boudart that the TMC's show a Pt-like behavior in sev- 
eral catalytic reactions^ According to more recent inves- 
tigations, "early TMC's and TMN's often demonstrate 
catalytic advantages over their parent metals in activ- 
ity, selectivity and resistance to poisoning" and "for sev- 
eral reactions, such as hydrogenation reactions, catalytic 
activities of TMC's and TMN's are approaching or sur- 
passing those of group VIII noble metals" m The TMX 
surfaces are also technologically important as substrate 
materials in growth processes, e.g., in wear-resistant mul- 
tilayer coatings on industrial cutting tools" and for 
growth of carbidic nanostructuresi^ 

For catalytic applications, the stable or ideal surfaces 
are not necessarily the most suitable ones. Often the 
best site for a reaction is found on a stepped or in some 
way non-perfect surface, e.g., at kinks or around de- 



fects, where less stable faces of the material are exposed. 
Such sites often host surface states or surface resonances 
(SR's). This calls for studies on surfaces that present 
such surface states or resonances. 

In previous studies, the reactivities of the TiC(lll) 
and TiN(lll) surfaces are attributed to the presence 
of SR's of both Ti and C/N character ^i . Calculated 
trends in atomic adsorption strength are explained with 
a concerted-coupling model (CCM), in which the atomic 
frontier orbitals interact with both types of SR's. More 
recently, we have indentified a descriptor, defined as 
the mean energy of the TM-derived SR (TMSR), for 
atomic and molecular adsorption and for activation- 
energy barriers Ai Hence, the existence of several linear 
relations between the atomic and molecular adsorption 
energies, "scaling relation" , have been shown^- Such re- 
lations are of importance in the design of novel types of 
catalysts^ - — 

This paper is devoted to a deeper and more generalized 
understanding of the chemisorption on the TMC(lll) 
surfaces, thus also laying the ground for the above men- 
tioned descriptor. We extend the work done in our pre- 
vious study of the atomic adsorption on TiC(lll) and 
TiN(lll) and use density-functional theory (DFT) to in- 



vestigate whether the CCM is applicable to other TMC's. 
Our method consists of a detailed study of the trends in 
reactivity along periods and groups of the substrate par- 
ent metal and of the adsorbate, correlated with a careful 
mapping and analysis of the underlying details of the 
changes in surface electronic structure upon adsorption. 

The substrates chosen in our study are the TMC's 
formed with the parent metals Sc, Ti, V, Zr, Nb, Mo, 
Ta, and W in NaCl structure (see Fig.Q]). As adsorbates 
we choose atomic H, B, C, N, O, and F as well as the 
molecules NH, NH2 and NH3. These particular choices 
of TMC's and adsorbates allow us to capture changes 
in adsorption properties along both periods and groups 
in the periodic table. All studied TMC's adopt a NaCl 
structure either in stable or in metastable phase. 



II. TRENDS IN SURFACE PROPERTIES 

In this Section the computational details and results 
from our DFT calculations on the clean TM-terminatcd 
TMC(lll) surfaces are presented. First, the stability of 
the (111) surfaces is considered by comparing the cleav- 
age energies of these surfaces with the corresponding re- 
sults for the (100) surfaces. Then, the relaxed surface 
structures are presented and compared with existing re- 
sults. This is followed by a charge transfer analysis and a 
detailed analysis of both the energy- and space-resolved 
surface densities of states (DOS's). Particular empha- 
sis is put on the presence and character of SR's on the 
TMC(lll) surfaces. Throughout the presentation, anal- 
yses of the trends with respect to the TM component of 
the studied TMC's are made. 
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FIG. 1: The early transition-metal carbides under investiga- 
tion. 



Experiments show that after heating and ion or elec- 
tron bombardment under low-temperature and low- 
pressure conditions the (111) surfaces of TiCji 6 - VCp 7 - 
ZrC^ NbC^^ and TaC^ 2 - are unreconstructed and 
TM terminated. Therefore, our study deals with un- 
reconstructed and TM-terminated (111) TMC surfaces. 
Surface electronic structure characterizations with angle- 
resolved photoemission studies (ARPES) show the pres- 
ence of surface states on the (111) surfaces of TiCj22r— 
ZrC, 26 NbC, 27 and TaC. 28 Also, experiments show 
that both H 2 and 2 dissociate on TiC (11 1)^29-32 
ZrC(lll)^r^ and NbC(lll)^ 6 -^ Despite this large 
number of experimental studies, theoretical investiga- 
tions on TMC(lll) surfaces are scarce4°-~— 

This paper presents our results in the following or- 
der. First, a trend study of the clean TM-terminated 
TMC(lll) surfaces is presented in Section [TTJ In par- 
ticular we address the properties of the surface elec- 
tronic structures that are relevant for the adsorption 
mechanisms. Then, in Section IlIII the adsorption-energy 
trends are presented together with the trends in the 
adorption-induced changes in surface electronic struc- 
ture. The different trends in adsorption-energy and 
electronic-structure changes are analyzed and discussed 
in Section [Vj thus connecting our results to the previ- 
ously proposed CCM. The paper is concluded by Sec- 
tion IVI1 which summarizes the main conclusions of our 
investigation and makes some outlooks to possible rami- 
fications and further investigations. 



A. Computational details 

The surface calculations presented in this paper 
are performed within the DFT formalism using the 
plane-wave pseudopotential code Dacapo. 44 The ion- 
electron interaction is treated with Vanderbilt ultra- 
soft pseudopotentials4£ The exchange-correlation energy 
is included by the generalized gradient approximation 
(GGA) using the PW91 functional^ 6 . We utilize a slab 
geometry, with slabs of 4 to 8 bilayers (a bilayer being a 
unit of one TM layer and one C layer) , a vacuum region 
thickness corresponding to 5 bilayers, that is, at least 
10.8 A, and periodic boundary conditions. Each atomic 
layer is composed of one atom in a (1 x 1) geometry. The 
atoms in the three (four) outermost atomic layers on the 
TM-terminated side of a 4 (5 — 8) bilayer thick slab are 
allowed to relax until the sum of the remaining forces on 
all relaxed atoms is less than 0.05 eV/A, while the re- 
maining atomic layers are fixed at the bulk geometry. A 
Monkhorst-Pack sampling 4 ^ of 8 x 8 x 1 special fc-points 
and a plane-wave energy cutoff of 400 eV are used. The 
slab used to model the (111) surface is asymmetric, which 
gives rise to a discontinuity in the electrostatic potential 
at the cell boundary. This is corrected for by using the 
scheme in Ref. |48l 

To characterize the surfaces we utilize several elec- 
tronic structure tools. A Bader analysis is used to calcu- 
late the charge localization around individual atomsj 49 i 50 
Total and local, that is, atom-projected, DOS's for the 
surface bilayer are obtained by projecting the Kohn- 
Sham wave functions onto individual atomic orbitals and 
plotted as a function of energy (relative to the Fermi level 
Ep). To identify the surface specific properties that arise 
upon creation of the surface, the surface DOS's are com- 
pared to the bulk DOS's by studying the differences be- 
tween the two quantities. Information about the spatial 
localization of the surface-localized states is provided by 
analyzing the space-resolved surface DOS, that is, the 
Kohn-Sham wave functions. 



B. Cleavage energies 

The cleavage energy .E c ieav, that is, the energy needed 
to create two surfaces upon cleavage of a bulk structure 
along a specific crystallographic plane, is calculated as 



-Ecloav — -Eslab (n) ~ ni?bul] 



(1) 



where -E s i a b in) is the total energy of a slab with n TMC 
bilayers that exposes the two surfaces under investigation 
and 



EhuXk — Eslabin) — -E s l a b(n — 1) 



(2) 



is the bulk energy of one TMC bilayer, if n is sufficiently 
large. Thus, the cleavage energy corresponds to the sum 
of the surface energies of the two surfaces obtained upon 
cleavage. In the case of calculations on stoichiometric 
TMC (111) slabs, the cleavage energy is equal to the sum 
of the surface energies of the TM-terminated surface and 
of the C-terminated surface. 

Our calculated E c \ cav values for the TMC(lll) sur- 
faces, after relaxation of only the TM-terminated side of 
the slabs, are given in Table fl] Along a period -E c i eav 
shows a maximum for group IV, thus showing the same 
variations as the ones found for the bulk cohesive energies 
in our previous study^i Down a group, the variations in 
£-cicav are small but discernible and do not show any ap- 
parent correspondence to the bulk cohesive energies of 
Ref. iH 

According to the calculated E c \ cav values, TiC(lll) is 
the surface that requires the most energy to create among 
the considered TMC surfaces. Nevertheless, this surface 
is routinely grown by chemical-vapor deposition (CVD) 
under high temperatures as wear-resistant coating on in- 
dustrial cutting tools^^ The calculated -E c icav value for 
the TiC(lll) surface agrees well with those of previous 
DFT calculations^^ 

Compared to the TMC(IOO) surfaces, the TMC(lll) 
surfaces are found to have a higher -E c ieav value (see Ta- 
ble U) and are thus less stable. Also, the variations in 
^cicav are larger for TMC(lll) than for TMC(IOO) sur- 
faces. Both these properties can be attributed to the 
polar nature of the (111) surfaced 



C. Surface geometry 

The calculated relaxations of the four topmost layers 
of the TMC(lll) surfaces are given in Table fl] A pro- 
nounced contraction (10-30%) of the first interlayer spac- 
ing is found on all the surfaces. At the same time, the 
second interlayer spacing increases compared to the bulk 
separation. Further down into the slab the relaxations 
become smaller and the structure converges to the bulk 
structure after three bilayers. 

Most of the TMC (111) surfaces exhibit an alter- 
nating positive-negative relaxation similar to the one 



found in metalsi 5 - 5 - Compared to the close-packed par- 
ent metal surfaces, where the outer-layer relaxations are 
only a few percents of the corresponding bulk interlayer 
spacingS) 55 ' 56 the relaxations of the TMC(lll) surfaces 
are larger and comparable to other ionic polar surfaces. 57 
This indicates that the bonding character in TMC's has 
a significant ionic contribution, as discussed in Ref. l5ll 
The relaxation of the (111) surfaces is hence quite dif- 
ferent from the smaller rumpled relaxation found on the 
non-polar (100) surfaces. 

The percentual relaxations of the (111) surfaces, rela- 
tive to the bulk interlayer spacings, decrease down each 
group. Along the periods, the smallest percentual relax- 
ation is found for group IV in period 3d and group V in 
period Ad. The larger percentual relaxations of 5-MoC 
and WC can be attributed to the NaCl structure being 
a metastable phase for these compounds. The largest 
structural changes are found on the ScC surface. The 
variations in percentual relaxations are directly corre- 
lated to the cohesive energies^ a small cohesive energy 
gives a large relaxation. 

For the contraction of the first interlayer distance there 
are some experimental data, which agree qualitatively 
with our calculated values and in most cases even quan- 
titatively, as shown in Table [fl For TiC(lll), tight- 
binding- and DFT— calculations agree well with our 
values. Several of our results are also in qualitative agree- 
ment with the theoretical ones in Ref. |43[ A deviation 
is found, however, for the NbC(lll) surface, where ex- 
perimental studies show a contraction of both the first 
and the second interlayer distance, by 15.5% and 4%, 
respectively^ The first observation agrees with our cal- 
culated value but the second one differs qualitatively 
from the contraction seen in our and other first-principles 
studiesi 41 i 42 



D. Bader charge transfer 

In the bulk TMC systems there is a charge trans- 
fer from the TM to the C atoms. 51 Table [IT] shows our 
calculated Bader charge values for the TM-terminated 
TMC(lll) surfaces. Both the values relative to the free 
atoms (in units of electronic charge |e|) and the values 
relative to the bulk values (in percentages) are given. 

Compared to the bulk, an accumulation of charge oc- 
curs on the first surface bilayer, mainly on the first TM 
atomic layer but also in several cases on the topmost 
C atomic layer. An exception is ScC(lll), where the 
C layer gains more electrons than the Sc layer. For all 
the investigated TMC's there is a total of 40-50% re- 
duction in the surface bilayer ionicity compared to the 
corresponding bulk ionicity^i For the NbC(lll) surface 
our ionicity value for the surface Nb atoms (+0.86e) is 
in good agreement with the calculated one of +0.90e by 
Zhang et al£^ 

To the right along each period (except for 8-M.oC), 
the charge accumulation on the first TM layer atoms in- 



TABLE I: Relaxed surface energetics and geometry of the considered TMC's: Calculated cleavage energies -E c i cav , in J/m 2 , for 
the TM-terminated (111) and for the (100) surfaces; Perpendicular equilibrium distances dij, in absolute values (A) and relative 
to the bulk values (in parentheses), between atomic layers at the TM-terminated (111) surfaces (index 1: topmost surface layer; 
index 2: second surface layer, etc.). The bulk interlayer distances A and the available experimental surface interlayer distances 
^i2 P are a l so gi ven - The calculated results are extracted from the 8-bilayer slabs and differ by less than 1% from the ones from 
the 4-bilayer slabs. 
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d cxp 
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III 


ScC 


6.06 


- 


0.967 


-28.5%) 


1.558 (15.2%) 


1.303 (-3.6%) 1.375 (1.7%) 1.352 


- 


3d 


IV 


TiC 


11.29 


3.29" 


1.019 


-18.5%) 


1.394 (11.5%) 


1.196 (-4.4%) 1.278 (2.3%) 1.251 


0.87 (-30%)' 




V 


vc 


9.14 


3.21" 


0.961 


-19.8%) 


1.293 (7.8%) 


1.176 (-2.0%) 1.204 (0.3%) 1.202 


1.09 (-10%) c 




IV 


ZrC 


10.09 


3.10 a 


1.117 


-17.8%) 


1.513 (11.5%) 


1.299 (-4.3%) 1.378 (1.5%) 1.357 


- 


Ad 


V 


NbC 


8.92 


2.96" 


1.094 


-15.6%) 


1.355 (4.5%) 


1.305 (0.6%) 1.323 (2.0%) 1.297 


1.09 (-15.5%) d 




VI 


5-MoC 


5.87 


1.86° 


0.998 


-22.4%) 


1.378 (7.2)% 


1.256 (-2.2%) 1.376 (7.1%) 1.285 


- 




V 


TaC 


9.24 


3.06" 


1.137 


-12.1%) 


1.319 (2.0%) 


1.318 (1.9%) 1.321 (2.2%) 1.293 


1.14 (-12%) e 


5d 


VI 


WC 


6.14 


- 


1.030 


-19.6%) 


1.319 (4.1%) 


1.254 (-1.1%) 1.354 (6.8%) 1.268 


- 



a The cleavage energies for the (100) surfaces are adapted from 
Ref. l53las -B c i cav = 2_E sur f, since a stoichiometric (100) surface slab 
has two equivalent sides. 

'Ref.lll 

c Ref.[Ij. 



d Ref. 
E Ref. 



TABLE II: The ionicity, that is, the amount of charge in 
units of |e| relative to the neutral atoms, of the atoms in the 
first (TM) and the second (C) surface layers (positive val- 
ues = donated electrons, negative values = gained electrons) 
obtained from a Bader analysis. The percentual changes in 
ionicity compared to the bulk^i are also given (positive values 
= more electrons). 
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Surface 


TM 


C 










rel. atom rel. bulk rel. atom rel. bulk 




III 


ScC 


+ 1.27 


+18% 


-1.95 


+27% 


3d 


IV 


TiC 


+1.08 


+28% 


-1.77 


+19% 
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VC 


+0.86 


+39% 


-1.47 


+4% 




IV 


ZrC 


+1.15 


+32% 


-1.93 


+14% 


Ad 


V 


NbC 


+0.86 


+48% 


-1.58 


-4% 




VI 


S-MoC 


+1.39 


+29% 


-2.29 


+16% 




V 


TaC 


+ 1.07 


+45% 


-1.93 


-1% 


5d 


VI 


WC 


+0.86 


+46% 


-1.53 


-4% 



lar and unstable^ To counteract this, a surface charge 
can be induced, creating a neutralizing electric field. In 
Ref. [53, it is shown that for the (111) surface of a crys- 
tal with the NaCl structure, this surface charge should 
be equal to 50% of the bulk ionicity, which is what we 
observe. 

The extra charge in the surface bilayer will affect its 
electronic structure and be of importance for its adsorp- 
tion characteristics, as described below in Sections III El 
and IIIII respectively. 



E. Electronic structure 

In this Section, we investigate in detail the electronic 
structure of the clean TMC(lll) surfaces to gain an 
understanding of which characteristics that are surface 
unique and potentially important for the adsorption. 

1. Bulk characteristics 



creases, while it decreases on the C layer atoms. The 
same trends are found down each group. 

The calculated charge accumulation at the TMC(lll) 
surfaces has its origin in the polar nature of these sur- 
face. As mentioned above, the TMC's are partially ionic 
materials. Therefore there is a macroscopic electric field 
caused by the non-zero perpendicular dipole moment of 
the TMC bilayer, which makes the (111) surface po- 



For convenience and to facilitate the discussion of the 
TMC(lll) surfaces, we here provide a very short sum- 
mary, based on Ref. l51l of some bulk characteristics of 
the TMC's. 

The bonding in bulk TMC's has contributions fronr 
iono-covalent TM-C bonds, from TM-TM bonds, and 
from C-C bonds. For all considered TMC's, the bulk 
DOS's and band structures consist of a low-lying valence 
band (LVB), dominated by C(2s) states, an upper va- 
lence band (UVB), with contributions from C(2p) and 



TM(d) states, and a conduction band (CB) of mainly 
TM(d) character. The covalent TM(d)-C(2p) bonding 
states are found in the UVB. The main contribution to 
the UVB (CB) comes from C(2p) [TM(d)] states, thus 
indicating the partially ionic character of the bond. The 
UVB and the CB, positioned on each side of Ep, are 
connected by a non vanishing continuum of TM-TM 
and TM-C states (for ScC also C-C states) states. The 
C(2p)-C(2p) bonds are found in the lower part of the 
UVB. 

Towards the right along each period, both the UVB 
and the CB are shifted down in energy relative to Ep and 
the UVB becomes less C localized. Down each group, the 
bands are shifted towards lower energies and the UVB 
becomes more C localized. 



the DOS-difference plots (Fig. [3]): (i) one region of mixed 
TM and C character, positioned in the upper part of 
the UVB and just below the TMSR (between the shaded 
yellow and red regions in Fig. [3]); (ii) several peaks of 
mixed TM and C character, positioned in between the 
CSR peaks in the lower part of the UVB; and (iii) one 
region of mainly TM character, positioned in the lower 
part of the CB. These regions of negative peaks corre- 
spond to strong quenchings of the bulk UVB and CB 
peaks, respectively. 



3. Trends in surface characteristics 



2. Common surface characteristics 

The total and atom projected DOS's of the first TMC 
surface bilayer are presented in Fig. [5] A common fea- 
ture of the bulk and surface DOS's is that both have an 
LVB, a UVB, and a CB. The bulk DOS is recovered in 
the third surface bilayer. Despite these similarities, all 
surface DOS's differ sharply from their respective bulk 
DOS's. To easier identify the surface properties we study 
the difference between the surface DOS (Fig. [5]) and the 
bulk DOS (adapted from our bulk study^i), as illustrated 
for VC in Fig. |3] 

For all considered TMC (111) surfaces, surface specific 
features similar to those found on TiX(lll)iiSi^ can be 
identified, in particular: (i) a TM-localized surface reso- 
nance (TMSR) , positioned in the bulk pseudogap at or in 
the vicinity of Ep (shaded red area in Figs. [2] and |3]) and 
(ii) a more strongly C-localized UVB than in the bulk, 
with several C-localized surface resonances (CSR's) in 
the lower part of the UVB (shaded yellow area in Figs. [H 
and[3]). The energetical extent of the surface UVB is, 
however, largely similar to that of the bulk UVB. 

Both kinds of surface resonances consist of states local- 
ized at the surface that overlap energetically with bulk 
states, hence the term surface resonance. They can be 
identified as positive peaks in the DOS-difference plots 
and from real-space analyses of the Kohn-Sham wave 
functions (Fig. EJ). These show that the TMSR's cor- 
respond to unsaturated TM bonds that extend into the 
vacuum [see Fig. [3(c)] and point towards the fee sites 
[see Fig. [3Jd)] where the C atoms would be present if the 
bulk stacking in the (111) direction had been continued. 
Hence the TMSR are dangling bonds with a three-fold 
symmetry. The C-localized states in the surface bilayer 
rearrange to form new bonds; however, they do not be- 
come fully saturated, which results in C dangling bonds 
striving in between the first-layer TM atoms towards the 
vacuum [see Fig. |3(a-b)]. Coupling of these states to the 
second bilayer bulk states gives them the status of surface 
resonances. 

In addition, several negative peaks can be identified in 



The analysis of the difference between surface and bulk 
DOS's shows that as the group number of the surface TM 
constituent increases along each period, the position of 
the TMSR is shifted to lower energies relative to Ep, as 
expected from the filling of the TM d states. The TMSR 
filling increases, varying from a non-filled TMSR above 
Ep for group III, to a partly filled TMSR at Ep for group 
IV, to a filled TMSR just below Ep for group V, and a 
filled low lying TMSR for group VI. At the same time, the 
TMSR amplitude decreases, while the width is constant 
or increases slightly. 

A shift to lower energies and a decreasing amplitude 
are observed also for the CSR's, starting from ScC(lll), 
whose CSR's lie just below Ep. In addition, the ampli- 
tudes of the negative DOS-difference peaks in the upper 
part of the UVB and in the lower part of the CB decrease. 
For ScC, both the DOS-difference plots and the Kohn- 
Sham wave functions show more pronounced CSR's than 
on the other TMC's. 

The DOS-difference plots show that as the period num- 
ber of the surface TM constituent increases down each 
group, the positions of the TMSR and of the CSR's are 
more or less unaffected, while the TMSR amplitude de- 
creases and the width increases slightly. 



4- Connection with experiments 

Several of the TMC(lll) surfaces have been studied 
by means of angle resolved photoemission spectroscopy 
(ARPES), revealing surface- localized states, there de- 
noted "surface states" ,22 r 25 i 29 We associate these with 
our TMSR's, termed "resonances" because of their loca- 
tion in the pseudogap where the bulk DOS is non vanish- 
ing. The ARPES study, with an experimental resolution 
of 0.2 eV, places the position of these TMSR's in the T- 
point at -0.2 eV for TiC^^^ at -0.2 eV for ZrC^ at 
-0.7 eV for NbC, 27 and at -0.7 eV for TaC, 28 in perfect 
agreement with our results. 
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FIG. 2: Total and atom projected densities of states for the first surface bilayer of the considered TMC(lll) surfaces. The 
shaded yellow areas indicate the bulk UVB regions and the vertical dashed lines mark the positions of the CSR's. The TMSR's 
are marked with shaded red areas. 




FIG. 3: Difference between the surface and bulk DOS's for 
the first surface bilayer of VC(lll). The shaded areas and 
vertical dashed lines represent the same quantities as in Fig. [2] 
Also shown are representative real-space Kohn-Sham wave 
functions illustrating the existence of (a-b) CSR's localized on 
the first bilayer C atoms at —4.9 V and —3.7 eV, respectively 
and (c-d) a TMSR localized on the first bilayer V atoms. 
Figures (a-c) are side views and (d) is a top view of the surface. 
Larger gray balls correspond to V atoms and smaller black 
balls are C atoms. 



F. Origin of the SR's 



moving to the right along a period (see Section IIID[) is 
due to the filling of the TM d states and therefore of the 
TMSR. At the same time the amount of extra charge on 
the C layer, compared to the bulk, decreases. 

For the ScC surface the whole DOS is positioned at 
higher energies compared to the other TMC's. We recall 
that in bulk ScC, the C-C bonds are more pronounced 
than in the other TMC's.— In addition, the energy sepa- 
ration between UVB and CB is largest for this TMC and 
the Sc-localized CB is very high up in energy relative 
to Ep. Therefore the ScSR states are not populated. 
However, the surface still strives to reduce its polarity, 
which is done by a charge accumulation on the C atoms, 
resulting in very pronounced CSR's. 



III. TRENDS IN ATOMIC ADSORPTION 

So far we have discussed the properties associated with 
clean TMC surfaces. In this Section the focus is on the 
atomic adsorption on the TM-terminated TMC(lll) sur- 
faces. Atomic adsorption is one of the first necessary 
processes in reactions at surfaces. We perform two series 
of trend studies: one with respect to the adsorbate and 
one with respect to the substrate. As adsorbates we con- 
sider period 1 and period 2 atoms H, B, C, N, O, and F. 
The substrate trend involves the change of the TM atom 
in the TMC. In our case we use the carbides ScC, TiC, 
VC, ZrC, NbC, 5-MoC, TaC and WC, thus spanning four 
groups and three periods in the periodic table (see Fig.rjJ). 
The consequences of changing the non-metal atom have 
been investigated in previous studies, where adsorption 
on TiC(lll) and TiN(lll) was conducted.-^ 



The appearance of the TMSR's is a result of the break- 
age of the iono-covalent TM-C bonds that cross the (111) 
cleavage plane, which causes (i) anti-bonding TM states 
in the empty bulk CB to collapse into more TM atomic- 
like states, which lie at a lower energy than the CB and 
are positioned at or in the vicinity of Ep , and (ii) bond- 
ing C states in the upper part of the UVB to vanish as C 
atoms are removed from the surface. This is seen in our 
DOS-difference plots (illustrated by Fig. [3]) as negative 
peaks in the upper part of the bulk UVB region and in 
the lower part of the bulk CB region. A similar finding 
was reported for TiC (lll)£i2 The existence of TMSR's 
supports the charge accumulation picture obtained from 
the Bader analysis. 

In an analogous way, the changes in the lower part 
of the UVB DOS can be interpreted to be due to the 
breakage of the bulk C-C bonds upon formation of the 
surface, which results in the formation of CSR's. The 
increase of Bader charge on the C atoms in the surface 
bilayer supports the existence of CSR's and arises partly 
due to the influence of the extra charge on the TM surface 
bilayer on the C atoms. 

The increase in amount of extra charge on the TM sur- 
face layer, compared to the bulk, that takes place when 



A. Computational details 

The systems are modeled by four bilayers of TMC with 
a 3 by 3 atom geometry in the surface plane. We use a 
cutoff energy of 400 eV and a Monkhorst-Pack sampling 
of 4 x 4 x 1 k points. The high symmetry adsorption 
sites fee, hep, top, and bridge are considered. Both the 
adatom and the three top-most surface bilayers are al- 
lowed to relax in all directions. 

The adsorption energies -E a ds are defined as 



E* A ds 



■(E t 



slab+adatom 



-t-'frcc 



ad ale 



x), (3) 



where E[ rcc a datom is the energy of an isolated spin- 
polarized atom. It is known that the -E a ds values are 
sensitive to the choice of GGA flavor. To get an un- 
derstanding of the flavor dependence we have performed 
calculations on TMC(lll) systems with adsorbed H, N, 
and O atoms using the RPBE GGA functional^ The 
RPBE functional has been shown to give better results 
than the PW91 GGA functional for adsorption energies 
on TM surfaces,— whereas for bulk structure determina- 
tion in some cases it gives worse results than the PW91 



TABLE III: Calculated atomic adsorption energies in eV. The 
preferred adsorption site is fee, except for C, N, and O on 5- 
MoC and N and O on WC, which favor hep site. 



Period Group Surface (site) 


H B C N F 


III ScC (fee) 


2.97 2.92 4.52 5.36 7.68 6.97 


3d IV TiC (fee) 


3.60 5.68 7.87 7.86 8.75 6.92 


V VC (fee) 


3.29 5.83 7.54 6.85 7.31 5.79 



IV ZrC (fee) 3.62 5.56 7.73 7.79 8.73 7.00 

4d V NbC (fee) 3.45 5.89 7.61 6.92 7.41 5.85 

VI S-MoC (fee) 3.08 5.70 6.76 5.61 6.33 5.10 

S-MoC (hep) 2.81 5.60 7.11 6.19 6.54 5.09 

3.48 5.95 7.73 7.10 7.53 5.71 

■ ></ VI WC (fee) 2.95 5.84 7.23 6.12 6.40 4.95 

2.69 5.73 7.22 6.32 6.59 4.89 



V 


TaC (fee) 


VI 


WC (fee) 




WC (hep) 



TABLE IV: Calculated perpendicular distances in A between 
the adatoms and the TMC(lll) substrates for the systems 
given in Table [TTT1 



Period Group 


Surface (site) H B C 


N 


O F 


III 


ScC (fee) 1.25 1.55 1.25 


1.07 


1.09 1.30 


3d IV 


TiC (fee) 1.03 1.32 1.10 


1.02 


1.07 1.28 


V 


VC (fee) 0.99 0.97 1.11 


1.07 


1.14 1.38 



4d 



IV ZrC (fee) 

V NbC (fee) 

VI S-MoC (fee) 



1.01 1.35 1.11 1.03 1.08 1.31 
1.03 1.33 1.15 1.10 1.17 1.42 
1.16 1.36 1.25 1.30 1.46 1.62 



S-MoC (hep) 1.24 1.35 1.26 1.24 1.34 1.59 





V 


TaC (fee) 


1.00 1.32 1.14 1.08 1.14 1.36 


0(1 


VI 


WC (fee) 


1.06 1.33 1.17 1.14 1.21 1.50 






WC (hep) 


1.21 1.34 1.19 1.17 1.25 1.56 



functional^ According to our calculations on H, N, and 
O adsorbed on TMC(lll) surfaces, the RPBE functional 
gives consistently lower -E a d s values than the PW91 func- 
tional. A similar type of variation was found for ad- 
sorption on pure TM surfaces^ As the decrease in E^ds 
values for each adsorbate is the same on all the consid- 
ered TMC(lll) surfaces (0.32 eV, 0.45 eV, and 0.12 cV 
for the H, N, and O adsorbates, respectively), the E a d s 
trends with respect to the substrate are not affected by 
the choice of GGA functional. Also, as these changes are 
smaller than 0.5 eV, the changes in -E a d s trends with re- 
spect to the adsorbate are not significant enough to affect 
the key features of these trends. 



strate. The smallest variations in -E" a ds values are found 
on S-MoC. 

For each adsorbate, the trends in E a ds along the sub- 
strate periods [Fig. 0Jb)] show that as group IV — > VI, 
(i) the adsorption strength decreases for the adsorbates 
N, O, and F; (ii) the -E a ds variations are small for H and 
C; and (hi) the adsorption strength increases slightly for 
B. The strongest adsorption is found on group IV TMC's 
(except for B). 

On ScC, the adsorption trends show somewhat differ- 
ent behaviors than on the other substrates. Also, the 
variations in -E^ds within each substrate group are much 
smaller than the ones within each substrate period. 



B. Atomic adsorption energies 

To elucidate our results in a clear way we gather the 
calculated -E a d s values in Table |TTT] and Fig. S) The 
most favorable adsorption site is fee, except for some ad- 
sorbates on S-MoC and WC, which favor the hep site, 
closely followed by the foe site. The relaxed perpendic- 
ular distances between each adatom and the TMC(lll) 
substrates are given in Table IIV1 In the following, the 
discussion will focus on the fee site. There are three ar- 
guments for this choice. First of all the shape of the -E a ds 
trends does not change if the hep preference is consid- 
ered. Secondly, as we are interested in focusing on the 
changes in surface electronic structure due to electronic 
factors we want to reduce the effects arising from geo- 
metric factors. Thirdly, in Ref. it has been shown for 
TiC(lll) that a similar adsorption mechanism applies for 
adsorption in both fee and hep sites. 

When moving from left to right along the adatom pe- 
riod 2 [B->F, see Fig.ija)] we find (i) "M" -shaped E ads 
trends for all TMC's; (ii) grouping, that is, similar -E a ds 
values within each substrate group; and (iii) strongest 
adsorption either for O or for C, depending on the sub- 



C. Bader analysis 

The Bader analysis gives a fractional electronic charge 
transfer from the substrate to the adsorbate (see Ta- 
ble |Vj). This indicates that the adsorbate-TMC bond 
is partially ionic. Most of the charge that the adsorbate 
gains is from the nearest three surface TM atoms, except 
on ScC where the C atoms are the main electron donors. 

The variations in the charge transfer down the TMC 
groups are much smaller than those along the TMC pe- 
riods. This indicates a more similar adsorbate-substrate 
interaction within substrate groups than within substrate 
periods. The largest amounts of charge transfer are found 
for the adsorbates N and C. However, the most filled 
outer electron shell is found for F, which is the most 
electronegative of the adsorbates. In our previous work 
on TiX(lll) we found, based on Bader analysis and the 
adsorbate-induced changes in electron density, that F has 
the strongest ionic adsorbate-substrate bondi^&iii 

Overall, there is no clear correlation between the calcu- 
lated -E a ds and charge transfer values, implying that the 
charge transfer alone cannot explain the -E a ds trends. In 
the same way as in our previous studies on TiX(lll) ) ' 8 ' 10 ' 
we show in the following that there is a significant cova- 
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FIG. 4: Calculated adsorption energies -E" a ds for period 1 and period 2 adatoms in fee site on the considered TMC(lll) surfaces 
as a function of (a) the adsorbate atom and (b) the group number of the substrate TM constituent. 



TABLE V: Calculated Bader charge transfers, in units of 
electron charge |e|, from the TMC(lll) surface to the atomic 
adsorbate for the systems given in Table [Hll 



Period Group 


Surface (site) H B 


C 


N 


O F 


III 


ScC (fee) 0.74 1.05 


1.45 


1.58 


1.33 0.84 


3d IV 


TiC (fee) 0.64 1.09 


1.35 


1.35 


1.17 0.80 


V 


VC (fee) 0.58 0.76 


1.16 


1.21 


1.09 0.76 



id 



IV ZrC (fee) 

V NbC (fee) 

VI 5-MoC (fee) 



0.67 1.19 1.49 1.50 
0.59 0.92 1.29 1.33 
0.65 1.16 1.55 1.54 



5-MoC (hep) 0.58 1.10 1.51 1.55 



1.25 0.82 
1.18 0.80 

1.26 0.82 

1.27 0.82 



V TaC (fee) 0.63 1.04 1.44 1.45 1.25 0.81 

5d VI WC (fee) 0.70 1.58 2.14 1.96 1.50 0.88 

WC (hep) 0.61 1.42 2.10 1.97 1.48 0.85 



lent contribution to the adsorbate-TMC(lll) bond. 



D. Density of states 

To learn about, in particular, the covalent parts of the 
adsorption bond, this Section is devoted to the trends in 
the adsorption-induced electronic structure. More pre- 
cisely, the difference in the surface DOS before and af- 
ter adsorption (ADOS), that is, the adsorbate-induced 
changes in DOS are investigated. With this very useful 
tool we monitor trends with respect to both substrate 
and adsorbate. 



1. Common characteristics in ADOS for atomic adsorbates 

Figures [5] and |6] show the ADOS's for C and O, re- 
spectively, adsorbed on the considered TMC(lll) sur- 



faces. Figure [7] shows the ADOS's for the different ad- 
sorbates on the VC(lll) surface. The general form of all 
the ADOS's consists of negative peaks of exclusively TM 
d character at the location of the clean-surface TMSR's 
(the shaded red regions) and negative peaks (or minima) 
of C character at the location of the clean-surface CSR's 
(vertical lines in the shaded yellow regions). Positive 
ADOS peaks are observed below, in between, and above 
the various negative peaks. A more detailed analysis of 
the atom-projected ADOS's shows that all the TMSR 
states associated with the three TM atoms close to the 
fee site, where the adatom is adsorbed, are depleted. 



2. Trends in ADOS with respect to the substrate 

The variations in ADOS between the different sub- 
strates are illustrated by addressing the differences be- 
tween two representative example: C and O atoms, re- 
spectively. 

The ADOS for adsorbed C, Fig. [5l shows the pres- 
ence on all substrates of a pronounced positive adsorbatc- 
projected DOS (represented by the blue line) peak that 
is pinned just below the depleted TMSR region and that 
extends with a low-amplitude tail throughout the up- 
per part of the UVB. This peak overlaps with substrate- 
TM states and can therefore be identified as a bonding 
adsorbate-TMSR level. The corresponding antibonding 
level, of mainly TM character, is found above the de- 
pleted TMSR region. The ADOS in the UVB consists 
of a number of subpeaks, of mixed adatom- C, substrate- 
TM, and substrate-C character, that are located in be- 
tween the depleted CSR levels. Such a structure indicates 
that the bonding adsorbate-TMSR level has interacted 
with the CSR's in the UVB. 

The ScC substrate forms an exception to the above 
pattern. Here, the interaction with the CSR's in the 
UVB appears to be much stronger. Also, the overlap of 
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FIG. 5: Total and atom-projected ADOS's for the adsorbate C atom in fee site on the considered TMC(lll)surfaces. The 
shaded areas and the dashed vertical lines represent the same quantities as in Fig. [2] 
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FIG. 6: Total and atom-projected ADOS's for the adsorbate O atom in fee site on the considered TMC(lll) surfaces. The 
shaded areas and the dashed vertical lines represent the same quantities as in Fig. [2] 
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the adsorbate-projected DOS peak with the substrate- 
TM states is not as pronounced as on the other TMC 
substrates. 

The ADOS for adsorbed O, Fig. El shows a much more 
varied structure. On several substrates, the adsorbate- 
projected DOS is now located inside the UVB region and 
shows a stronger substrate- C character than in the case 
of adsorbed C. When moving from left to right along each 
TMC period, this peak shifts to higher energies relative 
to the substrate UVB region: for group IV TMC's (TiC 
and ZrC) it lies at the lower edge of the UVB, whereas for 
group VI TMC's (S-MoC and WC) it lies at the upper- 
edge of the UVB. At the same time, its form changes: 
quite localized with a tail towards higher energies for 
TiC and ZrC; broad for VC, NbC, and TaC; and local- 
ized with a tail towards lower energies for <5-MoC and 
WC. These observations indicate a stronger and much 
more varied interaction between the bonding adsorbate- 
TMSR level and the substrate CSR's than in the case of 
adsorbed C. 

Again, ScC stands out from the above pattern, be- 
ing characterized by a much stronger participation of the 
substrate CSR's in the bonding. 

It can also be noted that for both adsorbed C and 
adsorbed O, the key features of the ADOS's are the same 
within each substrate TM group. 

These results show the presence of significant covalent 
bondings between the adsorbate and the substrate SR's. 
In particular, both types of SR's (TMSR's and CSR's) 
appear to participate in the chemisorption. This indi- 
cates that the previously reported picture for adsorption 
on TiC(lll) and TiN(lll), that is, a concerted-coupling 
model (CCM) should be valid on these TMC(lll) sur- 
faces as well. Such a picture is further pursued in Sec- 
tion [V] where the differences and variations in electronic 
structure are related to the trends in calculated -E a ds val- 
ues described in Section UlIBI 



3. Trends in ADOS with respect to the adsorbate 

To illustrate the variations in ADOS upon change of 
adsorbate, Fig. [7] shows the calculated ADOS's for the 
period 2 adatoms on the VC(lll) surface. This trend 
has previously been studied on the TiC and TiN (111) 
surfaces and described with the CCMi^r— 

For all the adsorbates there are quenchings of both 
the TMSR and the CSR's, but to different degrees for 
different adsorbates. As we go from left to right along 
the adatom period (B— >F), the energy of the adatom- 
projected DOS decreases. At the same time, its width 
changes: localized for B and C, delocalized for N and 
O, and again localized for F. Also, the depletion of the 
TMSR decreases gradually. The degree of depletion of 
the substrate CSR's varies also between the different ad- 
sorbates. An analysis of the Kohn-Sham wave functions 
shows that adatom-projected DOS peaks that lie in the 
lower part of the UVB consist of strong adatom-C bond- 




-4 -2 

E-E p (eV) 



FIG. 7: Total ADOS for the adsorbates B, C, N, O, and F 
on the VC(lll) surface. The shaded areas and the dashed 
vertical lines represent the same quantities as in Fig. [2] 



TABLE VI: Results from the calculations on the molecules 
NIU (x = 1,2,3) adsorbed on the VC(lll) surface: adsorp- 
tion energies -E a ds, perpendicular distances d between the 
molecule and the surface, and charge transfers from the sur- 
face to the molecule obtained by a Bader analysis. 





S ads (eV) 


d{k) 


Bader 


(units of 


e|) 


NH 


5.93 


1.23 




0.96 




NH 2 


4.63 


1.48 




0.64 




NH 3 


0.84 


1.83 




0.14 





ing states, while peaks in the upper part of the UVB 
contain strong adatom-TM bonding states. 

Again, these variations indicate strongly varying de- 
grees of interaction between adsorbate, TMSR, and CSR 
levels. They resemble the trends previously found on the 
TiC and TiN (111) surfaces and should therefore be pos- 
sible to explain in a similar way within the CCM. This is 
done in Section [V] where the DOS variations are related 
to the calculated E a d s trends. 



IV. MOLECULAR ADSORPTION 

In this Section we present results for molecular ad- 
sorption on the VC(lll) surface, which is here chosen as 
a prototype for molecular adsorption on the TMC(lll) 
surface. The molecular adsorbates are the NH X molecules 
(x = 1,2,3), which are a part of our study in Ref. [Ill 
The same supercell size and computational parameters 
as those used for the atomic adsorption are employed for 
these calculations (see Section Hi Ap . The molecules are 
adsorbed with the N atom closest to the surface in the 
fee site (the stable site for N atomic adsorption). The 
adsorption energy is calculated relative to the energy of 
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FIG. 8: Total and atom-projected ADOS's and LDOS's for N, NH, NH 2 and NH 3 adsorbed on VC(lll). The shaded areas 
and the dashed vertical lines represent the same quantities as in Fig. [2] 



an isolated NH X molecule as 



E, 



ads 



-(E t 



slab+NH T 



,)• (4) 



The charge transfer between substrate and absorbed 
molecule is calculated as the sum of the Bader charges 
on each of the constituent atoms in the molecule. 

A summary of the results is given in Table IVI1 As the 
molecule size x increases, the perpendicular distance to 
the surface increases, while both the adsorption energy 
-E a ds and the charge transfer from surface to molecule 
decrease. 

Figure [S] shows calculated LDOS's (projected on the 
molecule atoms) and ADOS's (projected on the substrate 
atoms) for N and NH, (x = 1,2,3) on VC(lll). The 
energy and structure of the LDOS's vary between the 
different adsorbates. For atomic N, it consists of one 
delocalized region in the middle of the UVB. For NH, 
there is a slightly weaker delocalized region of N character 
in the middle of the UVB, together with a sharp peak 
of mixed N and H character in the lower part of the 
UVB. For NH2, the delocalized region of N character in 
the middle of the UVB is further weakened and a new 



sharp peak of mixed N and H character appears below the 
UVB. For NH3, the delocalized region in the middle of the 
UVB disappears, while the sharp peak below the UVB 
is stronger and lies at lower energy. For all adsorbates, 
there is a quenching of both the TMSR and the CSR's. 
These SR depletions are largest for N and decrease every 
time a H atom is added to the molecule. Hence, each 
addition of a H atom weakens the interaction between 
the adsorbed molecule and the SR's. 



V. DISCUSSION 

As described above, we find several trends in the cal- 
culated adsorption energies: "M" -shaped variations with 
respect to atomic adsorbate, linear variations with re- 
spect to substrate, qualitatively different types of adsorp- 
tion on ScC and for F adatom, and trends in molecular 
adsorption strength. 

In this Section, we describe how these trends can 
be described within our previously proposed concerted- 
coupling model (CCM) for atomic adsorption on 
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TiX(lll) surfaces,!*^ First, the CCM is described by 
giving an account of the main results of Rcfs. 7.9.10 and 
shown to apply also on the here considered TMC's by 
pointing out the similarities in electronic structure re- 
sults. Then, the different adsorption-energy trends are 
explained in terms of the CCM. Finally, the results lay 
the ground for the single descriptor £ccm for the adsorp- 
tion strength, recently introduced in Ref. I 111 . 



A. Concerted-coupling model 

For TM surfaces, the ci-band model yields a successful 
description of electronic structure and adsorptioni 61 i 62 
For instance, its key parameter Ed, the mean energy 
of the substrate d band, is a good descriptor for, e.g., 
adsorption^r— Such a fact facilitates the design of new 
TM catalysts by computational screening*^— 

For adsorption on TMC's, however, there are devia- 
tions from the d-band modeL 42 ' 66 Indeed, we find that 
different TMC surfaces can have the same value of Ed 
but different -E^ds values. For example, the Ed values for 
TiC(lOO) and TiC(lll) surfaces are as close as 0.43 eV 
and 0.30 eV, whereas the -E a ds values for oxygen differ 
considerably, being 5.79 eV and 8.76 eV, respectively^ 

In the d-band model, adsorption trends are explained 
by the interaction between the adatom frontier orbital 
and the narrow TM d band. Using the terminology of 
the Newns- Anderson (NA) model for chemisorption^2r— 
this bond is typically "strong" and results in the for- 
mation of separated bonding and antibonding adatom- 
substrate states (see case d in Fig.O^ 

In Refs.@@, and[l0|, atomic adsorption on the TiC and 
TiN (111) surfaces is described as a result of two types 
of interactions between adsorbate and substrate. In the 
terminology of the NA model, the coupling of the adatom 
frontier orbital is typically "strong" with the substrate 
TiSR and "weak" with the substrate XSR's (X = C or 
N) (as in case d and in cases a-c in Fig. [SJ respectively). 
A concerted action of these couplings gives a qualitative 
explanation for the calculated adsorption energy trends 
on the TiC(lll)i2iifl and TiN(lll)^ria surfaces. The 
TiSR is present on the TiX(lll) surfaces but not on the 
(001) ones. The XSR's are found in the TiX(lll) UVB's. 

In the first mentioned coupling above, the large over- 
lap of the localized TiSR with the adatom orbital causes 
a strong adatom-TiSR interaction (in the NA sense). 
Well-separated bonding and antibonding states of mixed 
adsorbate and Ti character are then formed (case d in 
Fig. [SJ. The bonding-state energy lies below the free- 
adatom and TiSR levels, while the antibonding state re- 
sides above the TiSR level. 

In the other coupling, the bonding adatom-TiSR level 
interacts with the XSR's present in the substrate UVB. 
Due to the short range of the XSR's (compared to the 
TiSR), this interaction is weak (in the NA sense) and 
causes a mixture of broadening and shifting of the bond- 
ing adlevel-TiSR state. A state located in the middle of 



the UVB is mainly broadened (case b in Fig. [9]), while a 
state at the edge of the UVB is mainly shifted away from 
the UVB center of mass (cases a and c in Fig. ^. 

Evidence for this concerted coupling is given by 
detailed analyses of the calculated adsorbate-induced 
DOS's (ADOS's) and of real-space visualizations of the 
Kohn-Sham wave functions. These show that upon ad- 
sorption, there are (i) a sharp decrease in DOS at the 
TiSR energy; (ii) a sharp increase of DOS just above 
the TiSR energy; (hi) depending on the adatom species, 
narrow or broad bands of mainly adatom character at 
the edge of or within the substrate UVB energy range, 
respectively; (iv) a depletion of X-localized UVB states 
at the XSR energies; and (v) a formation of adatom- 
localized sub-peaks in between the energies of the sub- 
strate XSR's. Points (i)-(iii) show that bonding and an- 
tibonding adlevel-TiSR states are formed, while points 
(iii)-(v) prove a coupling between the bonding adlevel- 
TiSR state and the XSR's in the substrate UVB. 

More detailed analyses of the trends in calculated 
ADOS's and Kohn-Sham wave functions for second- and 
third-period adatoms on TiX(lll) surfaces show that (i) 
the magnitude of the DOS reduction at the TiSR energy 
decreases successively as the adatom number Z increases 
along a period and (ii) the adatom-X bonding character 
of the adatom-localized peaks increases as the adatom 
number Z increases along a period. These trends, which 
arise from the successive lowering of the adlevel energy 
as Z increases, indicate that the adatom-TiSR coupling 
decreases in strength along each adatom period, towards 
the right in the periodic table, while the contribution of 
the XSR's to the bonding increases. 

Such trends provide a basis for understanding the cal- 
culated trends in atomic adsorption energies. For in- 
stance, the maximum for the group- VI adatoms (O and 
S) is explained in terms of a stronger coupling to the sub- 
strate CSR's. On the other hand, the weaker bonding of 
group- VII adatoms (F and CI) is explained by a weak- 
ened coupling to the CSR's, due to the adatom state be- 
ing almost fully ionized by its interaction with the TiSR. 
Thus, the F and CI adsorbates lack almost any covalent 
interaction with the substrate UVB and their adsorption 
is practically ionic in nature, as confirmed by Bader and 
charge-density analyses. Also, the almost equal adsorp- 
tion strengths of C and N can be explained in terms of 
the opposite trends in adatom-TiSR and adatom-CSR 
coupling strengths found for varying Z values. 

On the same basis, the CCM is also able to describe 
changes in adsorption energies arising from changes in 
adsorption site, i.e., between fee, hep, and top O adatoms 
onTiC(lll)^ 



B. Generalization of the CCM to other TMC's 

The calculated energetics and electronic structures of 
atomic adsorption on the TMC(lll) surfaces described in 
Section IPlTl can be analyzed in the same way as described 
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FIG. 9: A schematic representation of the solutions of the Newns-Anderson model for atomic chemisorption on a metal 
surface.— The left and right panels illustrate the "weak" and "strong" chemisorption limits, respectively. The solutions, that 
is, the adatom-localized DOS's after adsorption, are yielded by the intersections between the straight lines a-d and the function 
A(e). The energy of the adatom frontier orbital before adsorption is given by the intersection of the lines a-d and the a;-axes. 
Thus, cases a and c correspond to a shift of the adlevel to lower and higher energies, respectively, while case d corresponds to 
a formation of well-separated bonding and antibonding levels. In case b, only a broadenin g of the adlevel is obtained. The 
function A(e) represents the DOS of the clean surface before adsorption. Adapted from Ref. |68| . 



above for TiX(lll). 

The fact that upon adsorption the TMSR's and CSR's 
are quenched on all considered surfaces shows that both 
types of resonance participate in the substrate-adatom 
bond. The strong character of the adatom-TMSR bonds 
is evidenced by the formation of well-separated bonding 
and antibonding states of mixed adatom and TM char- 
acter, at energies below and above the TMSR's, respec- 
tively. The weak character of the adatom-CSR bonds is 
confirmed by the broadening and/or shift of the bonding 
adatom-TMSR level. 

The exception is ScC(lll), on which the adatom-CSR 
coupling approaches a strong NA character, with a strong 
depletion of CSR states in the middle of the UVB and the 
formation of separated bonding and antibonding states 
below and above the UVB, respectively. This is due 
to the stronger CSR's on this surface, which in turn 
arises from the stronger C-C bonds in bulk ScC than 
in the other TMC's£± Also, the adatom-TMSR coupling 
is qualitatively different on ScC(lll), compared to the 
other TMC's, as the TMSR of clean ScC(lll) lies above 
Ep and is thus empty before adsorption. 

Thus, the basic adsorption mechanism is the same on 
all the considered TMC(lll) surfaces, implying that the 
general applicability of the CCM includes TMC(lll) sur- 
faces. In addition, ScC(lll) provides a playground for 
testing the contribution of the adlevel-CSR coupling to 
the total adsorption strength. 

In the following, we give more explicit support to the 
applicability of the CCM to the TMC(lll) surfaces by 
showing that the CCM can be used to understand the 



major features of the calculated atomic and molecular 
adsorption-energy trends on the considered TMC(lll) 
surfaces. 



Adsorption trends with respect to substrate 
and CCM 



Figure fjjb) shows the calculated energy values for 
atomic adsorption on TMC(lll) surfaces as functions of 
the substrate TM group number. In this subsection we 
show that the CCM can capture its main trends. 

In a simple molecular two-level picture, the strength of 
a bond is related to the energy shift between the levels 
of the isolated constituents and of the hybridized states 
resulting from the coupling. In a simplified version of 
the Newns model, an adsorbate and the active part of 
the substrate can be viewed as such a "molecule" ££■ Ac- 
cordingly, from the calculated ADOS's the bond strength 
can be extracted as the energy difference between the fi- 
nal bonding adatom level (i.e., after coupling to both 
TMSR and CSR's) and the TMSR peak in the surface 
DOS. The CCM is then able to explain the trends in bond 
strength by relating this energy difference to concerted 
adatom-TMSR and adatom-CSR interactions. 

In the following, we employ such an approach to un- 
derstand the qualitatively different trends in E^ds values 
for C and O adatoms, as the substrate TM group number 
increases from IV to VI [Fig. BJb)]. 

For C, the -E a ds values are approximately constant be- 
tween groups IV and VI. The calculated ADOS's (see 
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Fig. [5]) show that the final bonding adatom level lies 
at the upper edge of the UVB on all considered sub- 
strates. Thus on all these substrates, the energy differ- 
ence between the final adatom level and the original sub- 
strate TMSR is approximately constant. This implies 
that the result of the concerted action of the adatom- 
TMSR and adatom-CSR's interactions is approximately 
constant. In particular, Fig.[S]shows that the DOS's of all 
substrates resemble the same NA type of adlevel-CSR's 
interaction, that is, type c in Fig. [5J 

For O, with its lower 2p level, the E a a s values decrease 
strongly as the substrate TM group number increases 
from IV to VI. This is reflected in the calculated ADOS's 
(Fig. [6]) by a strongly decreasing energy difference be- 
tween the final adatom level and the original substrate 
TMSR peak. Thus, in contrast to C, for O the concerted 
action of the two types of coupling yields qualitatively 
different results. In particular, Fig. [5] shows DOS's that 
correspond to adlevel-CSR's interactions that vary in a 
qualitative way from (cf. Fig. [9J NA case a (on group 
IV), to NA case b (on group V), and to NA case c (on 
group VI). As a consequence, the adlevel-CSR's inter- 
actions cause a downward energy shift (on group IV), a 
broadening (on group V), and an upward energy shift (on 
group VI) of the bonding adlevel. 

We also notice that for both adsorbates the antibond- 
ing adatom-TMSR level lies above Ep (and is thus 
empty) for all substrates except the group- VI TMC's 
(5-MoC and WC. This filling of antibonding states im- 
plies an extra decrease in adatom-TMC bond strength, 
as group V — > VI. Indeed, this can be seen in the calcu- 
lated E ads trend for C [Fig. gfb)] . 

For all adsorbates except F, the calculated E^s val- 
ues increase from ScC(lll) to TiC(lll). As described 
above, the adlevel-CSR's interaction is more pronounced 
on ScC(lll) than on the other TMC's, thus approaching 
the strong NA limit. Also, the adatom-TMSR coupling 
is qualitatively different on ScC(lll), due to the TMSR 
of clean ScC(lll) lying above Ep and being thus empty 
before adsorption. The weaker adsorption on ScC(lll) 
can thus be understood to be due to the qualitatively 
different interaction with the two types of SR's. 



D. Adsorption trends with respect to adsorbate 
and CCM 



The adsorption-energy trends with respect to adsor- 
bate [Fig. Q|a)] show "M "-shaped E ad s trends for most 
of the considered TMC's. The trend on TiC(lll) was an- 
alyzed in Ref. and explained within the CCM to arise 
from competing opposite trends in adsorbate-TMSR and 
adsorbate-CSR's interactions. 

The calculated ADOS's for the here considered TMC's 
(illustrated for VC in Fig. [7]) show trends that are sim- 
ilar to those found for TiC. At the same time, the -E a ds 
trends for all considered TMC's resemble the one on 
TiC (111), with a local maximum for the O adatom. How- 



ever, on group-V TMC's the calculated E ad s values for 
O are approximately equal to the ones for C, and on 
group- VI TMC's the O values are smaller than the ones 
for C. Thus, the chemisorption strength of O appears to 
weaken relative to C as the substrate TM group number 
increases. This can be understood to be a consequence 
of the different E ad s trends for C and O that were de- 
scribed in Sec. IV CI above, with the adsorption strength 
decreasing much faster for O than for C as the TM group 
number increases from IV to VI. 

On the other hand, the .E a ds values for ScC(lll) show 
a monotonic increase between B and O adatoms. As de- 
scribed above, the adlevel-CSR's interaction is stronger 
on ScC(lll) than on the other considered TMC(lll) sur- 
faces. Therefore, the monotonically increasing trend (as 
B — > O) of the adlevel-CSR's coupling contribution to 
the adatom-TMC bond is stronger on ScC(lll). 

Finally, on all substrates there is a sharp decrease in 
adsorption energy, as O — > F. This arises from the weak- 
ened adlevel-CSR's interaction, which is due to the fact 
that the adatom is almost fully ionic from its interaction 
with the TMSR, as described previously on TiC (11 1)7 



E. Molecular adsorption and CCM 

The adsorption energies of the considered molecular 
adsorbates NH^ (x — 1, 2, 3) on VC(lll) show a de- 
creasing trend as a; increases (Table IVI[) . This trend can 
be understood in terms of the CCM as arising from the 
concerted coupling between each molecular orbital and 
the two types of SR's. 

The quenchings of both TMSR and CSR's in the cal- 
culated ADOS's for NFL, on VC(lll) (Fig. [5]) indicate 
the presence of couplings to both types of SR's. In ad- 
dition, the adsorbate-projected LDOS's show strongly 
bound sharp peaks of mixed N and H character as well as 
more delocalized regions of only N character within the 
UVB. A comparison with calculated DOS's for the free 
NH X molecules shows that the sharp peaks correspond to 
low-energy molecular levels that interact weakly with the 
substrate SR's. On the other hand, the delocalized re- 
gions arise from the coupling of higher-energy molecular 
levels with the SR's, in a way that is similar to the ad- 
sorption of atomic N. Figure [5] shows that as the number 
x of H atoms in the molecule increases, the interaction of 
the higher-energy molecular levels with the SR's weakens, 
thus lowering the adsorption strength. 



F. Applicability of the CCM 

The above discussion shows that for the atomic and 
molecular adsorption on the TMC(lll) surfaces a picture 
based on the concerted action of two types of adatom- 
substrate interactions applies. It can be used to describe 
key features of the calculated adsorption-strength and 
electronic-structure trends. In this concerted-coupling 
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FIG. 10: Correlation between _E a ds and eccM for the atomic 



model (CCM), two types of SR's participate in the bond, 
TM-localized SR's (TMSR) and C-localized SR's (CSR). 
Therefore we conclude that the CCM is valid for atomic 
adsorption on all the considered TMC(lll) surfaces. 



G. Descriptor for adsorption on TMC(lll): eccM 

To be able to describe the variations in adsorption 
strength in a simple yet efficient way, a descriptor 
£ccm > defined as the mean energy of the TMSR, can be 
introduced^ Due to the approximately constant energy 
difference between TMSR and CSR's in the considered 
compounds, such a descriptor is able to capture the im- 
portant variations of both TMSR and CSR energy trends 
while being at the same time both conceptually simple 
and measurable or calculable. 

This is confirmed by Fig. 1101 which shows a linear cor- 
relation between our calculated -E a ds values and £ccm 
for each of the considered atomic adsorbates, also when 
the results of Ref. [Tl| are significantly extended. The ex- 
ception is ScC(lll), which is a consequence of the men- 
tioned qualitative difference in electronic structure of this 
surface, that is, an empty TMSR and a strong adlevel- 
CSR coupling. As discussed above, the variations in the 
gradients of the lines for the different adsorbates can be 
understood within the CCM from the details in the in- 
teractions between the adlevel and the different SR's. 

In Ref. [ll|, linear correlations are also found between 
£ccm and the molecular adsorption energies of NH, NH 2 , 
and NH3 , as well as the activation-energy barriers for N2 . 
This is a consequence of the fact that both atomic and 
molecular adsorptions appear to follow the same basic 
mechanisms of the CCM, as argued above. As a conse- 
quence, we show in Ref. [II] that scaling and Br0nsted- 
Evans-Polanyi (BEP) relations apply for adsorption on 
TMC(lll) surfaces. Such findings are of importance for 



the design of new catalysts. 



VI. CONCLUSIONS 

The possibility to understand materials has today in- 
creased considerably, as DFT has developed into a predic- 
tive theory. One obvious application of DFT is to calcu- 
late numbers for bonding strengths, structure parameters 
and coordinates, etc. Another, slightly more demand- 
ing, application is to look for more detailed information 
and develop conceptual frameworks and models in which 
we can formulate our understanding and on which we 
can base our further thinking, including ideas about new 
materials. We have found that various densities of states 
(DOS's) are excellent tools for such understanding and 
vehicles for thought. In our approach, extensive use of 
total, local, projected, and difference DOS's are corner- 
stones in the building-up of a conceptual framework and 
model. 

This study deals with understanding the atomic and 
molecular adsorption on the TMC(lll) surfaces using 
such a detailed electronic structure analysis approach. 
By extensive DFT calculations on ScC, TiC, VC, ZrC, 
NbC, (5-MoC, TaC, and WC, we study trends in clean sur- 
face properties, adsorption energies, and various DOS's. 
In brief, we find that only a certain part of the surface- 
localized electronic structure, corresponding to the TM- 
localized (TMSR) and the C-localized (CSR) surface 
resonances (SR's) are of importance to understand the 
trends in adsorption energies from one TMC surface to 
another and from one adsorbate to another. Despite the 
TM termination of the investigated TMC surfaces, the 
second-layer C atoms are found to play a crucial role in 
the chemisorption. This is particularly evident on the 
ScC(lll) surface, where the TMSR is empty and where 
the CSR's are particularly strong. 

Having thus identified the key parameters, we possess 
the foundations for a concerted-coupling model (CCM) 
in which trends in adsorption strength are the result 
of a concerted action of both adsorbate-TMSR and 
adsorbate-CSR's couplings. This has earlier been shown 
for the TiC and TiN (111) surfaces!&±2 but the breadth 
and versatility are here shown by applications to other 
TMC's and trends, including: adsorption trends with 
respect to substrate, adsorption trends with respect to 
atomic adsorbate, adsorption trends on ScC(lll), ad- 
sorption trend for group- VII adatoms, and molecular- 
adsorption trends. This deepens the interpretation of 
the model and broadens its usefulness. 

The applicability of the CCM to both atomic and 
molecular adsorption opens up the possibility to study 
reactions for, e.g., catalytic applications on this class of 
materials. For example, it allows the formulation of a sin- 
gle descriptor eccm for the adsorption strength.— Also, it 
implies the existence of scaling relations between molecu- 
lar and atomic adsorption strengths as well as Br0nsted- 
Evans-Polanyi relations. 11 



18 



Since the CCM framework is based on rudimentary 
bonding principles with general applicability, we be- 
lieve that it is possible to generalize it to other mate- 
rials that possess surface-localized states. We have al- 
ready shown its applicability to atomic adsorption on 
TiN(lll)&iJi and believe that the same chemisorption 
mechanism should be valid for other nitrides. Ligand and 
vacancy systems have been shown to belong to the group 
of materials where the CCM applies, 11 as do certain TM 
surfaces*^ where, however, it does not need to replace 
the sufficient and natural d-band model. Natural exten- 
sions should include TM oxides, sulfides, and borides. 
Design of materials, including atomic-scale engineering, 



is also an enticing prospect for further applications. 
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